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A theory of Kondo lattices is applied to studying possible magnetic and charge structures of 
itinerant-electron antiferromagnets. Even helical spin structures can be stabilized when the nesting 
of the Fermi surface is not sharp and the superexchange interaction, which arises from the vir- 
tual exchange of pair excitations across the Mott-Hubbard gap, is mainly responsible for magnetic 
instability. Sinusoidal spin structures or spin density waves (SDW) are only stabilized when the 
nesting of the Fermi surface is sharp enough and a novel exchange interaction arising from that 
of pair excitations of quasi-particles is mainly responsible for magnetic instability. In particular, 
multiple SDW are stabilized when their incommensurate ordering wave-numbers ±Q are multiple; 
magnetizations of different ±Q components are orthogonal to each other in double and triple SDW 
when magnetic anisotropy is weak enough. Unless ±2Q are commensurate, charge density waves 
(CDW) with ±2Q coexist with SDW with ±Q. Because the quenching of magnetic moments by 
the Kondo effect depends on local numbers of electrons, the phase of CDW or electron densities 
is such that magnetic moments are large where the quenching is weak. It is proposed that the so 
called stipe order in cuprate-oxide high-temperature superconductors must be the coexisting state 
of double incommensurate SDW and CDW. 
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I. INTRODUCTION 

It is a long standing and important issue to explain 
two types of magnetism, local-moment magnetism and 
itinerant-electron magnetism, in a unified theoretical 
framework. Some of their physical properties are sim- 
ilar to each other and others are different from each 
other. Even if they are phenomenologically similar to 
each other, their microscopic explanations can be totally 
different from each other. For example, the spin suscep- 
tibility obeys the Curie- Weiss law in both types of mag- 
nets. The Curie- Weiss law in insulators is due to the tem- 
perature dependence of local thermal spin fluctuations. 
This mechanism is of leading order in 1/d, with d spa- 
tial dimensionality. On the other hand, the Curie- Weiss 
law in metals is a controversial issue. Two mechanisms 
were proposed: the temperature dependence of the.mode- 
mode coupling between intersite spirLfhictuationsEl^El and 
that of Weiss' magnetic mean fieldsotl The former is of 
higher order in 1/d, while the latter is of leading order in 
1/d. It is interesting which is responsible for the Curie- 
Weiss law in actual d = 3 dimensions. There is a sharp 
peak at the chemical potential in the density of states in 
typical ferromagnetic metals, and there is sharp nesting 
of the Fermi supJjace in typical antiferromagnetic metals. 
Recent theoriesErB showed that in such itinerant-electron 
magnets the mode-mode coupling plays a totally nega- 
tive role in the Curie- Weiss law. On the other hand, a 
novel exchange interaction arising from the virtual ex- 
change of pair excitations of quasi-particles can have a 
large temperature dependence consistent with the Curie- 
Weiss law;Cra Weiss' mean fields are given by magnetic 
polarizations multiplied by the exchange interaction. 

Magnetic structures are different between the two 
types of antiferromagnets. Helical structures are sta- 
bilized in typical local-moment magnets p and sinu- 



soidal structures or spin density waves (SDW) in typ- 
ical itinerant-electron magnets. When ordering wave- 
numbers are incommensurate, there are two or more than 
two equivalent wave-numbers depending on lattice sym- 
metry. An interesting issue is which are stabilized, single 
or multiple SDW, in what conditions. 

Itinerant-electron magnets lie in the vicinity of the 
Mott-Hubbard metairinsulator transition. A combined 
theory of Hubbard'sB and Gutzwiller'sO implies that in 
metallic phases the density of states is of a three-peak 
structure, Gutzwiller's quasi-particle band at the chem- 
ical potential between the lower and uppet, Hubbard 
bands lying far from the chemical potentiaLEJ It is un- 
questionable that the Mott-Hubbard splitting exists even 
in metallic phases in the vicinity of the transition. Such 
metals must show a crossover between the two types 
of magnetism as a function of temperatures. Denote 
the energy scale of local quantum spin fluctuations by 
ksTfc- They must behave as local-moment magnets at 
T ^> Tk because local thermal spin fluctuations are dom- 
inant, and they behave as itinerant-electron magnets at 
T <C Tk because magnetic moments are quenched by 
local quantum spin fluctuations. When physical pheno- 
mana relevant to the Mott-Hubbard transition are exam- 
ined, first of all, local spin fluctuations should be accu- 
rately taken into account. 

Local spin fluctuations are rigorously considered in 
any single-site approximation (SSA) that includes all the 
single-site terms. Such an SSA is reduced to deteonining 
and solving selfconsistently the Anderson model,EJ which 
is one of the effective Hamiltonians for the Kandn prob- 
lem. Because the Kondo problem was solved,E3~t3 many 
useful results are available in clarifying single-site prop- 
erties in lattice systems. The three-peak structure dis- 
cussed above can be easily understood by the SSA theory 
or the mapping to the Anderson model. The quenching 
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of magnetic moments by local quantum spin fluctuations 
in lattice systems is mapped to that in the Kondo prob- 
lem, so that Tk discussed above is nothing but the Kondo 
temperature. A perturbative treatment of intersite terms 
starting from an unperturbed state constructed in the SSA 
is nothing but a theory of Kondo lattices. Because the 
SSA is rigorous for Landau's normal Fermi-liquid states 
in infinite dimensions (d — > +oo)J13 it can also be formu- 
lated as a 1/d expansion theory. 

One of the purposes of this paper is to study which are 
stabilized in itinerant-electron magnets, single or multi- 
ple Q structures, in what conditions. This paper is orga- 
nized as follows: In Sec. [n], the theory of Kondo lattices 
is applied to deriving Landau's free energy. What are 
studied in Sec. [11 are magnetic and charge structures of 
ordered states below the Neel temperature. Discussion 
is given in Sec. 0, and conclusion is given in Sec. 0. 
Landau's free energy for the Anderson model is studied 
in Appendix ^[ 



II. FORMULATION 



A. Theory of Kondo lattices 



Consider the following single-band model: 



H = -E%4^ +E Ufl ^ fl n 



2 /-^i ' *J 



(2.1) 



with hia = d' icr d t(T and n t = + n tl ; s lv = 

hY^aftAa^diP ( U = X i V and Z )> witn tne Pauli 

matrix. When only the first and second terms are consid- 
ered, the model (2.1) is reduced to the Hubbard model, 
in which the wellknown superexchange interaction is in- 
volved. Actual systems must be described by a multi- 
band model. A supjerexchange interaction arising from 
multi-band effectsML3 is phenomenologically included 
by the third term. The fourth term is the long range 
Coulomb interaction: = e 2 /|Rj — Rj|, with R, 

lattice sites. No spin-orbit interaction or no magnetic 
anisotropy is included in this paper. 

Assume Landau's normal Fermi-liquid state. The 
single-particle selfenergy is divided into single-site and 
multi-site terms: S cr (ie n ,k) = S CT (i£ n ) + AS cr (ie Tl , k). 
The multi-site term, AE cr (ie„,k), is of higher order in 
1/d, and is ignored here. The single-site term, E^ («£„), is 
of leading order in 1 jd. In the presence of infinitesimally 
small external fields, Ti' ext = — J2ia [A/i'+|(T5/jBff']n iff , 
it is expanded in such a way that E CT (i£„) = Eo + 
(1 - 4> 7 )ie n + (1 - 4> s )\agiJLBH' + (1 - 4> c )Ap' + ■ ■ ■ for 
I £ra | ^ ^bTk and T <C TV, with Tk the Kondo tempera- 
ture discussed in Introduction. Because («£■„) is equal 
to that of the mapped Andeipepjnodel, it follows accord- 
ing to Yamada and Yosida'sL!-ali3 that (f> 7 = ((f) s + (f> c )/2; 
and 



are approximately given by Eqs. (All) .and 
( A12j ) in Appendix. The electron density is given bytEl 



n d = 



N 



(2.2) 



with TV the total number of unit cells. When the system 
is in the vicinity of the Mott-Hubbard transition, n<j is 
close to unity so that (f> s 3> 1 and (f> c <C 1. The single- 
particle Green function is given by 



G a (ie n ,k) = J- — 



Vl iS n - £(k) ' 

for |e„| -C UbTk and T <C Tk, with 



i 



E - y t i:j e y 3 



(2.3) 



(2.4) 



The quasi-particle density of states is given by p*(e) = 
(1/N) J2 k 8(e — £(k)). When contributions from local 
spin fluctuations are only considered, the specific heat 
at T <C Tk is given by C — 7T + ■ ■ ■, with 7 = 
(2/3)^1^(0). 

When irreducible two-point polarization functions 
in the wave-number representation are denoted by 
TTa-r (itoi , q) , those in spin and charge channels are 
given by n s (iui,q) = n aa (iuji,q) - 7r ._ .(icJ;, q) and 
w c (iui,q) = n aa (iLUi,q) + -K a -^(iuJu q)- Spin and charge 
susceptibilities of the model (2.1) are given bycll 



Xs(m,q) 



2ir s (iuji,q) 



1- [t/+ij(f)(q)] 7r s (^,q) 



and 



Xc(iuJi,q) = 



2Tr c (iuJi,q) 



[U + 2V( c )(q)] ^(iwi.q)' 



(2.5) 



(2.6) 



with jW(q) = Y,j&i) jg ,) e- i *< R «- B *> and V^(q) = 
Y^j(M) Vpi c 'e~ l,i '( Ri ~ R i\ Irreducible polarization func- 



tions are also divided into single-site ones, Tr s (iuii) and 
TT c (iuji), and multi-site ones, A-K s (iuii, q) and A7T c (iuJi, q): 
ir s (iuji,q) = K s (iu>i) + Air s (iu t , q) and 7r c (iw,,q) = 
TT c (iu)i) + Air c (iuii,q). Because the single-site -terms are 
equal to those of the mapped Anderson modeljlil its spin 
and charge susceptibilities are given by 



and 



Xs(iui) 



Xc(i^i) 



2tt s (iuji) 

1 - UTT s (iuJl) 



2TT c (iLOl) 
1 + UTT c (iuJl) 



(2.7) 



(2.8) 



SIB 



The Kondo temperature is defined by ks 
l/[xs(0)] T _ >OK - When the Fermi- liquid relationtj'LJ is 
made use of, it follows that IzbTk — (0 7 /<^ s )//9*(O) so 
that IcbTk is as large as the bandwidth of quasi-particles. 
In this paper, our study is restricted to the strongly cor- 
related regime defined by 



k B T K /U « 1. 



(2.9) 



The Kondo temperature should be calculated by deter- 
mining and solving selfconsistently the mapped Anderson 
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model. However, it is treated as a phenomenological pa- 
rameter in this p aper . 

Eqs. (2.5) and (2.6) are also written in such a way that 



and 



with 



and 



Xc(«w ; ,q) 



1 - 3 / s («cj;,q)x s («w;) 

1 + V(iu>i,q)Xc(my 

j c s T ( P )f s , 2ATT s (iuJi,q) 
I s {lW h q) = J w (q) + . r 



V(iu h q) = V^(q) 



2ir c (iu!i)iT c (iu)i, q) 



It follows from Eqs. (|2.5|) and (|2.7|) that 



and 



7r s (0,q) = -[l + O(k B T K /U)} 



U0) = jj[l + O(k B T K /U)} 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



Here, U ^> \J i fHq)\ is assumed. To leading order in 
k B T K /U, Eq. ( [2~12j ) becomes 



I s (iuJ h q) = j( p Hq) + 2U 2 Air s (iLU,,q) 



It follows from Eqs. (2.6) and (2.8) that 



[U+V c (q)} 7r c (0,q) = 



±{U + Vc(q)} Xc(0,q) 
-i[C/+K(q)]Xc(0,q) 



(2.16) 



(2.17) 



and Un c (0) = |U"Xc(0)/[l - |t/x c (0)]. It is likely that 
[C/+V c (q)]|7r c (0,q)| > 1 and U\n c (0)\ > 1 in actual mag- 
nets. Note also th at <fi c <C 1 for rid ~ 1. Then, the 
second term of Eq. ( 2.1 3| ) can be ignored. In this paper, 
Eq. (2.13) is approximated by V(iu>i,q) ~ V^(q). In 
d = 3 dimensions, 



a 



Aire 1 



G| 



N ^ 



Aire 



N |p- G| 2 



(2.18) 



with G wave-numbers of reciprocal lattice points, v 
the volume of a unit cell, and the summation over 
p restricted to the first Brillouin zone. Because 
lim q ^G V^(q)/k B TK — +oo, the charge susceptibility 
vanishes for commensurate q: 



lim Xc(iui,q) = 0. 



(2.19) 



The susceptibility given by Eq. J2.1C|) with Eq. ( |2.16| ) 
is consistent with a physical picture for Kondo lat- 
tices that local spin fluctuations at different sites inter- 
act with each other by intersite exchange interactions. 



Then, I s (iu>i , q) should be called an exchange interac- 
tion, which is respo nsible for magnetic instability. The 
second term of..Ea. (2.16) is mainly composed of three 
contributions.E3Ej One is the superexchange interaction, 
which arises from the virtual exchange of pair excitations 
across the Mott-Hubbard gap. The phenomenological 
term J^^q) should be determined in such a way that 
the sum of this superexchange interaction and J^ p > (q) 
should be equal to the superexchange interaction of ac- 
tual systems. Another is the novel exchange interaction 
arising from that of pair excitations of quasi-particles. 
Whether it is ferromagnetic or antiferromagnetic depends 
on the dispersion relation of quasi-particles. For example, 
it is ferromagnetic in the so called flat-band and band- 
edge models.EJa On the other hand, it is antiferromagnetic 
when the nesting of the Fermi surface is sharp or when 
the chemical potential is at the band center. The other 
is the mode-mode coupling term among intersite spin 
fluctuations. Each of the superexchange interaction and 
the novel exchange interaction can cause magnetic in- 
stability, while the mode-mode coupling suppresses mag- 
netic instability. In the following part of this paper, we 
assume that the sum of the two exchange interactions 
is antiferromagnetic and is so strong to cause magnetic 
instability; the Neel temperature TV is determined by 
1/ \xs{0, Q)] t =t n = or 1 - i [1.(0, Q)x s (0)] t=Tn = 0, 
with Q ordering wave-numbers. 

As discussed in Introduction, local spin fluctuation 
are quite different between the two temperature regions, 
T Tk and T <C Tk- Then, local- moment magnetism 
is characterized by TV 2> Tk while itinerant-electron 
magnetism by TV <C Tk- In this paper, our study is 
restricted to itinerant-electron magnetism or TV <c Tk- 



B. Landau's free energy 

Magnetizations appear below TV. In this subsection, 
Landau's free energy in the presence of external fields 



-g^B 



^(H^-Si) -^A/xA, (2.20) 



is studied. For this purpose, it is convenient to consider 
Hn=H + E ext +H L , (2.21) 
instead of TL + Ti. ex ti with 

i i 

including no operator and 
H L = -J2 X t } i^-m w ) -Y, X< i C) (^-^) ( 2 - 23 ) 



including Lagrange's multipliers, \ ff and \\ c> . The ther- 
modynamic potential f2 for Tin is given by 



k (c) 



-n/fc B T = Tr exp 



k B T 



, (2.24) 
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with u the chemical potential. Lagrange's multipliers 
should satisfy 

on 



dX 



(<0 



and 



with 



on 



ox 



2<s<„)n = 



= i%i - (ni)n = 0, 



(2.25) 



(2.26) 



)n standing for thermal averages for Tin- Then, 
magnetizations and charge densities are given by sets of 
mi and iii, {m} and {n}. It follows that 



on 

dm iv 



X^ 



'({m},{n}) - -gn B H.u 



and 



^ = A l (c) ({m},M)-A^ 



(2.27) 



(2.28) 



Here, it is not necessary to consider the dependence of 
f2 on Jmj and {re } through A-^ and A- c ^ because of 
Eqs. ( |2.25| ) and When Q is minimized MS il func- 

tion of {m} and {nj, it follows from Eqs. ( [2.27 ) and 

( |2.28|) that X[ S J = \gii B H w and x\ c) = Au;, so that Hn 
becomes H + H ex t- The set of Eqs. fl2.25| ) and ( ]2.26D is 
nothing but the selfconsistency condition to determine 
magnetizations and charge densities for H. + H. ext , and f2 
defined this way is Landau's free energy for Ji + H ex t- 
The normal Hartree term is given by Ao = — ( \ Un^ + 

J2j ^ij C ' n <i) i with n d the average electron density defined 
by Eq. (|]^).E3 All of its effects are i ncluded in the unper- 
turbed state considered in Sec. [I A. Anomalous Hartree 



and Fock terms, which are due to spin or charge polar- 
izations, are given by 



A 



(«) 



i . 



-Urn; 



i E 4^ 



for spin channels, and 



A 



< c ) = x(c) 



1. 



a/ - ( -rA„, 



for the charge channel, with 
An, = n» — 



■ E ^ c) ^ 

i(^) 



(2.29) 



(2.30) 



(2.31) 



Here, Lagrange's multipliers, A^, and AJ , are included, 
and thermal averages (si„)n and (nj) n are replaced by 
mi V /2 and rif because Eqs. (2.25) and (2.26) are eventu- 

_(s) —(c) 

ally satisfied. When A it/ and A 4 are perturbatively con- 
sider ed s tarting from the Fermi-liquid state considered in 
Sec. |[IA|, i t follo ws fr om the selfconsistency condition of 
Eqs. (j2.25|) and (|2.26|) that 



2 E ^ ' ) A ^ + 2 E 7r " 2 cihh-h) A^i A 



M -(c) 

J2 



3U2 



- 2 E^ 

ilia 



Ah3i,32)X n X J2U 



+2 E E E 7r "" 1 ^ ^ ' j ' 2 ' j ' 3 ) 



Jl"l J2^2 J3"3 



xA (5) A (s) A (s) + 



(2.32) 



and 

An,i = 2 } V c 2(i, j)X) 

3 

+2 E E n ^ ii' ia^iii^L + ' ' ' > ( 2 - 33 ) 

iii2 "i 



with Tx v i{i,j) static two-point polarization functions 
in spin channels, ir c 2(i,j) those in charge channels, 
Vct'i,^,^), ^vcv(h,i2,h) and 7r C!/ 2 (i^ i 2 , i 3 ) static 
three-point ones with two spin vertexes and a charge ver- 
tex, and n VlV2V3Ui (11,12, ia, i±) static four-point ones with 
four spin vertexes. 

Note that 7r s (0,q) = £V (i, j) e -^< Ri - R ^ for v = 
x, y and z, and 7r c (0,q) = £\ (*> j) e - lc >-( R '- R ^ . It 
follows from Eqs. ( 2.3 2| ) and ( |2.33 ) that in the wave- 
number representation 



m q „ I \ - , , T (s) T (c) 

- x ^-^N q ^; A(q ' qi ' q2) v, 



x^ = 



Q2 



1 



7r s (0,q) 



E E 

V1V2V3 qiq2q3 



7T„ 



,(q,qi,q2,q3) 



and 



A (c) 



with 



An r 



xA (s) A (s) A (s) 
qi^i q2^ q3"3 



1 



(2.34) 



X ^qii/i^q 2 ^i "I ' (2-35) 



An c 



— An,e 
AT ^ 



A< s ) = 1 VA (s) e- iq ' R * 
q^ _/y A " / 



q jy 1 



T 0) 1 v^T (s:) -iq-Ri 

= a$ + + ^ (p) (q)"v, 



= A q =)-i[/An q -W c )(q)An q , 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



(2.41) 
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;e -i[(qi-R-ii)+(q2-R.« 2 )+(q3-Ri 3 )] j (2.42) 



and 



(,qi,q2,q3,q4 



xe 



-i[(qi-Ri 1 ) + (q2-Ri 2 ) + (q3-Ri 3 ) + (q4-R i4 )] 



(2.43) 



Because of the translational symmetry, 
TIVi 1^2 ^3 (01,02,03) oc (5 

qi+q2+q3 

and 

(qi,q2,q3,q4) oc <y qi+q2+q3+q4 , 

with 

' 1, q = G 
0, q^G' 



5 q = 



(2.44) 
(2.45) 

(2.46) 



Eqs. |[3j) and ( |2~35| ) can be solved in an iterative 
manner. When |m qi/ | < 1 and |Ara q | < 1, Eqs. ( |2.34| ) 

and ( p35| ) give trivial relations, A q J = w ql ;/Xs(0,q) an d 
A q ^ = An q /xc(0, q). When these approximate relations 
are made use of, it follows that 



A (s) - 
A q i/ — 



27r s (0,q) 2 



Vm„ v 



and 



-( C ) = Aw q 
q 2vr c (0,q)- 



(2.47) 



(2.48) 



In Eq. ( [2.47 ), Eq. (|2.14 ) is al so m ade use of. When these 
are used in Eqs. ( 2.34[ ) and ( |2.35 ), it follows that 



A q s J = 7n ~, \ ~ 2 E .9(q>qi>q2) m q^ A? 

XsV u >OJ qiq2 

+ E b uv , lU2U3 (q,qi,q2,q3) 



wi^2^3 qiq2q3 

X 771 qiZyi 77i q2l/2 ?72 q3 ,y 3 -|- 



(2.49) 



and 

A ( C) = An q 
q Xc(0,q) 



Y f(q,qi,q2)m qi 



i^iqiq2 



with 



and 



.g(qi,q2,q3) = 



^ 2 ^ c (qi,q2,q3) 
2 2 ^ c (o, qi ) 



(2.50) 



(2.51) 



i(Ol, 02, 03,04) 



' 2 3 



1(01,02,03,04)- 

(2.52) 



Eqs. ( [2.27 ) and ( |2.28| ) can also be written in the mo- 
mentum repre senta tion, and L agrange's multipliers are 
given by Eqs. ( 2.4S| ) and ( |2.50 ). By integrating them, we 
obtain 

O = + flext + ^s 2 + ^c 2 + ^s 2 c + ^s 4 + ' ' ' , (2.53) 

with flo the term including no spin nor charge polariza- 
tions, 



a 



-N 



\dVB ^2 (H q • m_ q ) + Y A Mq Ar 



fi s 2 = 



2Xs(0,q) : 



2Xc(0,q)' 



(2.54) 



(2.55) 



(2.56) 



fi s 2 c = -iV^ 5 (qi,q 2 ,q 3 )(m qi -m q2 )An q3 , (2.57) 
and 



qiq2q3 



iV v- 

= -j }^ ^ 4 (qi.q2,q3,q4) 



qiq2q3q4 



(m qi - m q2 )(m q3 - m q J - (m qi • m q3 )(m q2 • m q J 



+ (m qi m q4 )(m q2 -m q3 
In Eq. ( 2.58| ), the following relations of 

^ 4 (qi,q2,q3,q4) = ^1(01,02,03,04 



(2.58) 



,(qi,02,q3,04) 



(2.59) 



are made use of; b UlV2l/3Ui (qi, q%, q3, q4) = for other 
combinations of v^s. 

When ordering wave-numbers are incommensurate, 
there can be several equivalent ones to them. Denote 
the ith pair by ±Qi, and assume that not only ±Qi but 
also ±2Qi and ±4Qi are incommensurate: ±2Qi ^ G or 
±4Q, ^ G. When SDW with ±Q, are stabilized, CDW 
with ±2Q,; appear because of t he co upling term between 
SDW and CDW given by Eq. (§J57j) or 

fi s 2 c = -^^ 5 (Q l ,Q l ,-2Q i ) 

i 

x [(m Ql ) 2 An_ 2Qi + (m_ Q! ) 2 An 2Qi ] . (2.60) 

Here and in the following part, the summation over 
is made over pairs or it is made in such a way that one 
of ±Qj is only considered for each i. When only SDW 
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or helical structures with Qi and CDW with 2Qi are 
considered, 



&ext = -N^2<^gfi B (H Qi -m_ Q! ) + (H_ Qi -m Qi 

i ^ 

+ [A A i 2Qi A?i_ 2 Q ! + A/x_2 Qi An 2 Q 4 \, (2.61) 



^ 2 = N Y1 



\m Qi 



x.(o,Q<)- 



4- Xc (0,2Q 4 ) 



(2.62) 



(2.63) 



and 



i 

+2B 2 |(m Q! -m Qi )| 2 } 
+ |^{2B 3 |m Qi Hm Qj | 2 

+25 4 [|(m Qi - m Q3 )| 2 + |(m Qi . m_ Q .)| 2 ] 
+J B 5 [|(m Q , • m Qj )| 2 - |(m Qj • m_ Qj )| 2 ]}, 

(2.64) 

with 

S 1 = & 1/ 4(Q,-Q,Q,-Q), (2.65) 

B 2 = ^(Q,Q,-Q,-Q), (2-66) 



#3 = &„4 (Q, -Q, Q', - Q') + by* (Q, -Q, - Q', Q') 
-5MQ, Q', -Q. -Q') - ~MQ, -Q'< -Q, Q')> 

(2.67) 



S 4 = \b„4 (Q, Q', -Q, -Q') + i MQ, -Q', -Q, QO. 

(2.68) 

and 

B b = 2MQ, -Q. Q': -Q') - 2 MQ, -Q, -Q', Q')- 

(2.69) 

Here, Q are Q' (Q ^ ±Q') are two different wave- 
numbers among considered Q-j's, and symmetrical rela- 
tions such as 6„4(qi,q 2 ,q 3 ,q 4 ) = b v i (q 4 , qi, q 2 , q 3 ) = 
6 ! /i(— qi, — q 2 , q 3 , — q 4 ) and so on, are made use of. 

Four-point polarization functions are also divided into 
single and multisite terms so that 

Wqi,q2,q3,q4) = V + A6 i ,4(q 1 ,q 2 ,q 3 ,q 4 ). (2.70) 



The single-site term is approximately given by Eq. ( A10] ): 
&„4 = 2/x s (0) = 2ksTK- When a diagram correspond- 
ing t o Fi g. 2(6) of the previous papeiu is considered and 
Eq. (2.3) is approximately used, the multi-site term is 
calculated so that for T <C Tr- 



A6 l/ 4(q 1 , q 2 , q 3 , q 4 ) — 5 qi+q2+q3+q4 



xk B T K 



^ N ie »~^( k ) «e«-^(k+qi) 
1 1 
ze„-^(k+qi-(-q 2 ) (k+qi+q 2 +q3) ' 

(2.71) 

Here, the Ward relation for the three-point vertex func- 
tion in spin channels is made use of, as it was in the 
previous paper .□ In general, 



|B 5 | «Min(B 1) 5 2) B 3) B4)- 



(2.72) 



The nesting of the Fermi surface is characterized by 
£(k) ~ — £(k+Q) in a wide k region in the vicinity of the 
Fermi surface. When the nesting of the Fermi surface is 
sharp, B\ is much lager than other Bi's so that 



B x > Max(S 2 ,B 3 ,S 4 ) > \B 5 



(2.73) 



III. INCOMMENSURATE SDW AND CDW 

A. Sinusoidal or helical structures 

In this subsection, it is examined which are stabilized, 
sinusoidal or helical structures. The second order term 
in magnetizations, Af2 s 2 , does not depend on magnetic 
structures. When magnetizations are small, magnetic 
structures are determined by the fourth order term AQ S 4 . 
When single Q structures are assumed, the fourth order 
term is written as 



N 
~2 



{B 1 [2|m Q | 4 -|(m Q .m Q )| 2 ] 

+25 2 |(m Q -m Q )| 2 }. (3.1) 



When magnetizations are written as rriQ = rriq + i m'q, 
with rriQ and rriq being real, it follows that 

(m Q ■ m Q ) = (my 2 - (m Q ) 2 + z2(m Q ■ m Q ). (3.2) 

It is easy to see that |(mq • mq)| 2 becomes the smallest 
for helical structures, which are characterized by |m'q| = 
| m'q | and m'q _L niq, such as 



iq-rnqj 



0. 



(3.3) 



On the other hand, it becomes the largest for sinusoidal 
structures, which are characterized by niq || m'q, such 



as 
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|(m Q .m Q )| 2 Hm Q | 4 . (3-4) 
Therefore, helical structures are stabilized when 

B x - 2B 2 < 0, (3.5) 
and sinusoidal structures are stabilized when 



B 1 - 2B 2 > 0. 



(3.6) 



When no multi-site terms are included, B\ = B 2 = 
B 3 = B4 = b v i and B^ — 0. When there is no 
sharp nest ing of the Fermi surface, therefore, it is likely 
that Eq. (3.5) is satisfied. The novel exchange inter- 
act ioxL-jCan never be strongly antifcrromagnctic in this 
case.t£tt3 However, antiferromagnetic instability can still 
occur when the superexchange interaction is strongly an- 
tiferromagnetic. Even helical structures are stabilized in 
such itinerant-electron antiferromagnets. When the nest- 
ing of the Fermi surface is so sharp for incommensurate 
Q that Eq. ( |3.6| ) might be satisfied, sinusoidal structures 
are stabilized as is expected. 



B. Single or multiple Q structures 

In this subsection, it is examined whether multiple Q 
structures are stabilized or not. First, examine helical 



structures under an assumption that Eq. (3J3) is satisfied. 
Because (rriQ • hiq) = for helical structures, Af2 s 4 is 
reduced to 



Si|m Q ,| 4 + |^{2 J B 3 |m Q ,| 2 |m Qj . 



+2B 4 [\(m Qi ■ m Qj )| 2 + |(m Q , ■ m_ Qj .)| 2 ] 
+B 5 [|(m Qi -m Qj )| 2 - |(m Qt .m_ Qj )| 2 ]}. 



(3.7) 



When Eq. ( p.5[ ) is satisfied, it is likely that Bi < 2B 3 , 
B\ < 2B4, and I-B5I is much smaller than other Bi. In 
such a case, multiple helical structures are never stabi- 
lized. When helical structures appear at temperatures 
much lower than Tk, they must be of single Q. 

Next, examine sinusoidal structures under an assump- 
tion that Eq. (|3.6|) is satisfied. In this case, we can put 



JO. 



m, and m q ; 



-iOj 



nij, with nij being real. 



The fourth order term is given by 



N x , 



2B 2 



m 4 



+7V]T[i? 3 |m l | 2 K| 2 



2£?4(rrv m.j 



(3.8) 



When the multiplicity of Q is two or three, fi s 4 becomes 
the smallest when 



(mi 



m, 



= 



(3.9) 



for any pair of i ^ j. When Eq. ( |3.6| ) is satisfied, it is 
likely that 



B x + 2B 2 > 2B 3 



(3.10) 



is also satisfied. Then, multiple sinusoidal structures are 
stabilized and their polarization vectors are orthogonal 
to each other. It follows that 



with 



NX 



m2 N +-B ± (A)m 4 
X S (0,Q) 2 ±K ' 



B ± (\)=B 1 + 2B 2 + 2(\~1)B 3 



(3.11) 



(3.12) 



Here, |mi| is simply denoted by m because all of |nij| 
are the same as each other, and A is the multiplicity of 
ordering wave-numbers, A = 1, 2 or 3. Below T < Xjv, 



1 



and 



n. 



XN 



2B ± (X)\ Xs (0,Q)\ 2 ' 



(3.13) 



(3.14) 



Because Eq. ( ^3.14 ) is an increasing function of A, a sym- 
metry broken ordered state is never stabilized as long as 
A < 3; if there are three equivalent Q's, for example, 
single or double Q structures are never stabilized. 

When the multiplicity of Q is four or larg er than 
four (A > 4), it is impossible that Eq. (3.£) is sat- 
isfied for any pair of and Qj. Only two possi- 
ble magnetic structures are examined in this paper; 
triple or quartet Q structures. Triple structures are 
symmetry broken states, where only three Q's among 
A wave-numbers are ordered and their magnetizations 
are orthogonal to each other. Their free energy is 
given by Eq. (3.14) with A = 3. In quartet struc- 
tures, four Q components are ordered and each pair of 
magnetizations make the same angle, cos _1 (— 1/3) = 
121.63- •• degree; for example, mi = m (0,0,1), m 2 = 
to(0, v/8/3,-1/3), m 3 = m(V6/3,-\/2/3,-l/3), and 
m 4 = m (—- s/6/3, -V2/3, -1/3). Then, it follows that 



n S 2 + n si = an 



1 



Xs(0,Q 4 ) 2 



-Bom 4 



with 



B Q = B 1 + 2B 2 + 6B 3 + -B 4 . 
When the free energy is minimized, 

il s 2 + f2 s 4 = — 



2B Q | Xs (0,Q)| 1 



below T N . When A/B Q > 3/B ± (3) or 
B x > 4B 4 - 2B 2 + 2B 3 , 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



is satisfied, quartet Q structures are stabilized. Other- 
wise, triple Q structures are stabilized. One can conclude 
that quartet Q structures can be stabilized only when the 
nesting of the Fermi surface is extraordinarily sharp. 
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C. Coexistence of SDW and CDW 



IV. DISCUSSION 



The coupling betwee n heli cal str uctu res and CDW van- 
ishes according to Eqs. ( 2.60 ) and (3J5), a nd th at bet wee n 
SDW and CDW exists according to Eqs. fl2.60| ) and 
Once SDW with ±Q appear, therefore, CDW with ±2Q 
appear. The amplitude of CDW is given by 



An_ 2Q = g(Q, Q, -2Q) Xc (0, 2Q)m 2 Q . 



(3.19) 



When multi-site terms are ignored, <?(Q, Q, — 2Q) is ap- 
proximately given by Eq. ([A9|) : 



<?(Q,Q,-2Q) 



A/tt, n d < 1 
-A/tt, n d > 1 



(3.20) 



with A the hybridization energy of the mapped Aj 
derson model. According to the mapping condition 
1/A ~ -Im(l/iV)^ k G ! CT (+iO,k), so that the coupling 
constant is rather large. However, the charge suscepti- 
bility is small. Then, the amplitude of CDW must be 
small. In particular, it va nishes when 2Q is commensu- 
rate because of Eq. (2.19). 



Denote polarizations as hiq = 5 me 3 and A? 



2Q 



I At 



= , with m being real. Then, the coupling term 



is given by 



Atl s 2 c = --Ng(Q, Q, -2Q)xc(0, Q) 
x|m| 2 |An|cos(26» s - 6 C ). 
The free energy takes its minimum when 



2Z?r, n d < 1 
[ (2/ + 1)tt, n d > 1 ' 

with I being integers. Then, it follows that 

nii = mcos(Q • + S ) 

and 

|An|cos[2(Q-R i +6' s )], n d < 1 



29 x — 6 r 



Am = 



-|An|cos[2(Q-Ri + s )], n d > 1 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



The quenching of magnetic moments by local quantum 
spin fluctuations, which is one of the most essential effects 
in Kondo lattices, sensitivel y de pends on local number of 
electrons. According to Eq. (A2) in Appendix, the Kondo 
temperature at the ith site is approximately given by 



2A 



K 



1 — (n<j + Arij)! 



(3.25) 



Eqs. ( 3.23| ), ( |3.24 ) and ( 3.25| ) show that electron numbers 
Arij are modulated in such a way that magnetic moments 
are much quenched where Tk are high. In the less than 
half- filled case (n d < 1), for example, doped holes go 
mainly into sites where magnetic moments are small. 



Because no magnetic anisotropy is taken into account, 
directions of magnetic polarizations are totally indepen- 
dent of those of wave-numbers Q in this paper. In ac- 
tual magnets, absolute directions of magnetizations are 
mainly determined by magnetic anisotropy. However, 
the relative angles of magnetizations between different 
Q components must be 90° or at least close to 90° in 
double or triple Q structures even if magnetic anisotropy 
is taken into account; the magnetic part of the free energy 
is lower when the relative angles are closer to 90°. 

In Kondo lattices, the coupling between SDW and 
CDW is strong and the phases of SDW and CDW are 
never independent of each other. Then, even nonmag- 
netic impurities have a large pinning effect on coexisting 
SDW and CDW. 

Incommensurate SDW were observed in many metal- 
lic magnets, some of which . are of hiffikxaiough symme- 
try such as cubic CeAl 2 ,E-aE3 cubic Cr,t3c3 cuprate-oxide 
high-temperature superconductors, c3 which are approxi- 
mately regarded as orthorhombic lattices, and so on. In 
this paper, Landau's free energy is obtained up to the 
fourth order in magnetizations. Because saturated mag- 
netic moments are large in CeAl2 and Cr, higher than 
the fourth order terms are required to discuss their phys- 
ical properties at T -C Tn. However, it is still interest- 
ing to examine whether or not a multiple Q structure is 
actually stabilized in CeAl2 because magnetic moments 
are small just below T/v. It should be mentioned tl 
triple Q structure was actually proposed for CeAl 2 J 
The transition at Tjv — 38°C in Cr is of first order, and 
it is not certain if magnetizations are small enough even 
just below_T/v. Because a multiple Q structure was also 
suggested,U it is interesting to reexamine SDW in cubic 
Cr. It has been claimed that the so called_stripe or- 
der must be stabilized in the cuprate oxidesEj Because 
magnetizations are small, the theory of this paper is ap- 
plicable to magnetic states in the cuprate oxides. Within 
the theoretical framework of this paper, the stripe order 
must be nothing but the coexisting state of incommen- 
surate SDW and CDW. It is interesting which is actu- 
ally responsible for two equivalent satellites in neutron 
diffraction, two equivalent magnetic domains whose vol- 
umes are accidentally almost the same as each other or 
a double Q structure of SDW and CDW. 

The spectral weight of Gutzwiller's quasi-particle band 
is small in strongly correlated electron systems in the 
vicinity of the Mott transition. A large part of the spec- 
tral weight exists in the lower and upper Hubbard bands, 
which are far from the chemical potential. The formation 
of SDW and CDW causes not only gaps in quasi-particle 
spectra but also pseudogaps in the lower and upper Hub- 
bard bands. Their pseudogaps are as large as f7|m(Q)|, 
although gaps in quasi-particle spectra are of the order of 
ksTpf. It is interesting to observe changes in the density 
of states caused by the formation of SDW and CDW not 
only in the vicinity of the chemical potential but also far 
from the chemical potential. 
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V. CONCLUSION 

Magnetic and charge structures in strongly correlated 
electron liquids are studied within the theoretical frame- 
work of Kondo lattices. When there is no sharp nest- 
ing in the Fermi surface and the superexchange interac- 
tion is strongly antiferromagnetic, helical structures are 
stabilized even in itinerant-electron magnets. When the 
nesting of the Fermi surface is sharp enough, sinusoidal 
structures are stabilized. When incommensurate order- 
ing wave-numbers are multiple, in particular, multiple 
sinusoidal structures are stabilized. Their magnetic po- 
larizations are orthogonal to each other in double or triple 
sinusoidal structures, when magnetic anisotropy is small 
enough. 

Because the quenching of magnetic moments by local 
quantum spin fluctuations sensitively depends on local 
numbers of electrons in Kondo lattices, incommensurate 
CDW is inevitably driven by incommensurate SDW. In 
hole or electron doped systems at the vicinity of the Mott 
transition, therefore, the phase of CDW is adjusted in 
such a way that the density of doped carriers, holes or 
electrons, is large at sites where magnetic moments are 
small. 
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APPENDIX A: NON-MAPPED 
SINGLE-IMPURITY ANDERSON MODEL 

Consider the Anderson model with infinitely large U, 
constant hybridization energy A, and conduction band- 
width 2D. According to a previous paperj£3 the ground 
state energy of such an Anderson model is given by 



Er 



e d -»-\c-J (fc B T K ) 2 e-VA + J- A 2 , (Al) 



with ksTx = V DA exp [Tr(e d — /Lt)/2A]. Here, e d is the 
energy level of strongly correlated electrons and /i is the 
chemical potential; A c and A„ (v = x, y and z) are exter- 
nal fields or Lagrange multipliers, which are determined 
in such a way that magnetic moments and electron num- 
bers are m = (m x , m y .m z ) and n. 

In the absence of fields, the electron number is given 



by n d = - [dE G /dX c ]- x 



A c =0,A„ 



1 - (ttT k /2A), so that 



kBT K 



2A 



(l-n d ). 



(A2) 



In the presence of fields, magnetic moments and the elec- 
tron number are given by 



dE G 



I£e-*A./A + E>/A 2 



(A3) 



and 



8Eg 
d\c 



1 



T 2 -ttAc/A 

i K e 



2A 



(A4) 



It follows from these equations that 

2A m„(l - n) 



A l/ (m, n) 



ami 



A c (m,n) 



2A 



■In- 



1 — m 2 



1 — n 



7T (1 - n d )Vl - m 2 
Define a thermodynamic potential by 

w(m, n) = i?c( m 7 n ) + ™A C + m u \ u 



(A5) 



(A6) 



2A 

e d - n (1 - n) 

7T 



1-n 



2A 

.(l_ n )f^l n 

7T (1 — nd)vl — m 2 



(A7) 



It is easy to confirm that d£j(m, n)ldm v = A„(m, n) and 
duj(m,n)/dn = A c (m,n). 

The case when the electron number is smaller than 
unity is examined so far. The case where the electron 
number is larger than unity can also be treated within 
this theoretical framework if the hole picture is taken. 

In the presence of actual external fields, Afi and H, the 
thermodynamic potential to be minimized as a function 
of m and n is given by 



w(m, n) - nAfi - ig/i B (m • H) 

e d - - T K - nA/i - ^g^ B {m ■ H) 



1 



2X,(0) 
1 

2Xc(0) 



m 



- m 

2 



An - gm An - 



(A8) 



with An = n- n d , x s (0) = (tt/2A)/|1 - n d \, Xc(0) = 
(7r/2A)|l-n d |, and 



a = 



A/tt, n d < 1 
-A/tt, n d > 1 



(A9) 



Here, the results for n d > 1 in the hole picture are inter- 
preted into those in the electron picture. Then, it follows 
that 

2 4A 

6 - 4 = —77^ = — I 1 - n d\ = 1k B T K . (A10) 

When the Friedel sum ruleU is made use of, the den- 
sity of states at the chemical potential is given by 
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p(0) = sin 2 (7rnd/2) /ttA. When the Fermi- liquid rela- 
tions, x s (0) = <j> s p(0) and Xc(0) = <ft c p(0), are made use 
of, it follows that 



(™ d /2) 2 



|l-n d | sin^(7rn d /2)' 



and 



7 ,, i (™d/2y 

4>c = |1 - n rf|-7-27 



sin (nnd/2) 



(All) 



(A12) 
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